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ABSTRACT 

A  new  proof  is  given  for  Hausdorf f  *  s  condition  on  a  set  of  orients  wh i 
determines  when  the  function  generating  these  moments  is  in  L^.  The  proof 
uses  Legendre  polynomials  and  their  discrete  extensions  found  hv  Tchehvc he  * . 
Then  an  extension  is  given  to  a  weighted  L“  space  usina  Jacobi  polynomials 
and  their  discrete  extensions. 
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SIGNIFICANCE  AND  EXPLANATION 


The  paper  describes  a  method  of  obtaining  in  terms  of  the  moments 
approximations  to  the  solutions  of  the  finite  moment  problem 

1 

(1)  /  f(x)xVdx  =  U  ,  (v=  0,1,2,...)  . 

0  Y 

In  his  paper  [2j  Hausdorff  gave  conditions  on  the  moments  for  the 

problems  (1)  to  have  a  solution  f(x)  which  is  squares  integrable.  However, 
our  approximations  are  constructed  in  terms  of  the  coefficients  of  the 

Legendre  series  expansion 

00 

f(x)  -  i  c^ Py(2x  -  1),  (0  <  x  <  1)  , 

v=0 

where  P^(x)  are  the  Legendre  polynomials.  The  main  result  is  that 

Hausdorff s  condition  for  a  square  integrable  f(x)  are  here  expressed  in 

terms  of  the  c  .  This  transition  from  the  u  to  the 
v  v 

a  set  of  orthogonal  polynomials  on  the  discrete  set  x 
due  to  Tchebychef. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


cy  is  done  by  usinq 
0,  x  —  1,...,X“  n 


HAUSDORFF'S  MOMENT  PROBLEM  AND  EXPANSIONS  IN  LFGENDRE  POLYNOMIALS 


k.  Askey,  I.  J.  Schoenberg,  and  A.  Sharma 

1,  Introduction.  We  refer  to  [3]  for  a  description  of  the  problem  of  °ellman, 
Kalaba  and  Lockett  [1]  of  obtaining  approximation  to  the  inverse  Laplace 
transform.  They  reduce  the  problem  to  the  solution  of  the  finite  moment  problo- 

1  v 

(1.1)  /  f(x)x  dx  =  u  ,  (v=0,1,...,n-1; 

0  V 

and  obtain  approximations  for  f(x)  by  applying  Gauss'  n-point  quadrature  formula 
to  the  integrals  (1)  and  use  numerical  approximations  to  the  inverse  of  the  matr-> 
of  the  system  so  obtained. 

In  [3]  it  is  shown  that  the  inverse  of  the  Gauss  matrix  is  not  needed,  f’-'-er 
approximations  to  f(x)  are  obtained  if  we  determine  the  polynomial 

n-1 

(1.2)  R  . (x)  =  l  c  P  (1  -  2x) 

n-1  5  v  v 

of  degree  n-1  which  is  the  least  square  approximation  to  f(x)  in  [0,1] 

having  moments  ,u  ,...,p  „ .  The  coefficients  c  in  (1.2)  are  qiven  bv 

0  1  n-1  v 

lower  triangular  transformation 

v 

M-3)  cy  =  ( 2v  +  1)  l  (-1)1(V  *  *)(%  ,  (v  -  0,1 . n  -  1,  . 

i=0 

The  numerical  problem  of  Bellman,  Kalaba  and  Lockett  is  thereby  solved.  , 

this  approach  shows  that  the  infinite  problem  in  [0,1] 
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1 

(1.4)  /  f(x)x  dx  =  y  ,  (v  =  0,1,2  ...  to  infinity) 

0 


might  be  attacked  in  terms  of  the  Legendre  series  expansion 


(1.5) 


f  ( x)  -  y  o  P  (1  -  2x)  . 
«  v  \> 
v=0 


Hausdorff  devoted  to  the  problem  (1.4)  his  famous  paper  [2]  in  which  he  showed  the 
following: 

A.  The  system 

r1  v 

(1.6)  J  x  di|>(x)  »  u  ,  (v  *  0, 1 ,2, . . . ) 

0 


has  a  non-decreasing  solution  \(i(x)  if  and  only  if 

Any  =  y  -  (")y  +  •••  +  (-1)ny  >  0 

m  m  v1'  m+1  m+n 

for  m,n  >  0. 

B.  The  system  (1.6)  has  a  solution  ty(x)  of  bounded  variations  in  [0,1]  if 
and  only  if 


n 

v 

L 

V=0 


f" 

'■v 


)l 


.  n-v 
A  U. 


0(1)  as  n  ♦  ®  . 


For  a  direct  derivation  of  Hausdorff s  conditions  for  A  and  B  see  [4]. 

C.  The  system  (1.6)  has  a  solution 


x 

^(x)  =  /  <p(x)dx 
0 


where  £(x)  e  LP(0,1)  with  1  <  p  < 


if  and  only  if 


<n+i)p-1  [  ((")UB-Vl]p-  0(i; 


A  particular  case  of  C  is  this  (p  =  2): 


The  moment  problem 


(1.7) 


/  f(x)xVdx  =  u  , 
0 


(v  —  0,1,2,.*.) 


has  a  solution  f(x)  e  IT  (0,1)  if  and  only  if 


(1.8) 


S_  -  (n  +  1)  I  (")  (An'VUv)  =0(1) 
'*  u=0 


While  Hausdorff's  results  A  and  R  are  most  apt,  it  seems  that  the  result  (1.8)  might 


be  profitably  deduced  from  the  expansion  (1.5).  Since  /2n  +  1  r  (1  -  2x)  are 
orthonormal,  we  derive  from  the  Riesz-Fisher  theorem  and 


f ( x )  ~  y  - - —  /2v  +  1  P  (1  -  2x ) 

**  . -  v 

0  /2 V  +  1 


the  following:  The  moment  problem  (1.7)  has  a  solution  f(x)  e  L2(0,1)  if  and  only 


(1.9) 


l  2-  =  0(1 )  for  all  n 

v=0 


Formula  (1.3)  can  he  inverted  and  assumes  the  form 


(1.10) 


1  ,r2v+1^  f2v+1 i  .  .  v,2v+1 

v  K  *  (  v-iK . .  (  o 


(v+1)l  v  j 


Substituting  (1.10)  in  (1.8)  leads  to 


-3- 


These  simple  expressions  were  a  surprise  and  suggested  that 


S 

n 


n 

\ 

L 


v=0 


a 

n,  v 


2 

c 

v 


with  ^  given  by  a  reasonably  simple  expression.  This  will  be  shown  ir.  the  n*> 
section 


2.  Hausdorff  Theorem  via  Orthogonal  Polynomials 

Using  the  notation  introduced  in  the  first  section,  we  have 


.n-v  r  v,«  .n-v,, 

A  u  =  J  x  (1  -  x)  f(x)dx 

V  0 


n  1 

=  1  c.  i  P,  (1  —  2x)xV(1  -  x)n  Vdx 

k=0  k  0  k 


n  k  (-k)  .  (k  +  1 )  1 

=  I  ck  [  -1-  1T--2  /  xV+3<1  -  x,-'- 

k=0  j=0  ‘  M3'  0 


(2.1) 


Hv  +  1  )  r ( n  -  v  +  1  )  v  -k,  k  +  1,  v+1  1 

T(n  +  2)  ,  Ck’3F2'  1  ,  n+2  : 

k=n 


The  shifted  factorial  (a)n  is  defined  by 


(a>n  = 


T(  n  +  a  ) 
r<a) 


Murphy's  formula  for  Lecrendre  polynomials  was  used 

%<*>  r 
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To  obtain  (2.  *5),  write  ([4),  (4.1.3))  as  an  identity  between  hyperaeomotric  r.erier 


that  is 


a 


x)  = 


2F1 


-k ,  a 
a+1  -k-c' 


x) 


and  integrate  with  respect  to  a  beta  distribution.  Take  a=k+  1,  b  =  x  +  1, 
c  =  n  +  2  and  d  =  1  in  (2.5)  to  get 


(2.6) 


v*> 


(n  +  1  -  k)k 

~Tn~2)V~ 


,-k,  k+1 ,  -x 
3F2*-  1,  -n  ' 


1)  • 


Using  (2.6)  above  gives 


S 

n 


1 

n  +  1 


n 

r 

l 


k=0 


2 


(n  +  1  -  k)J 


l  KM*)]  • 

x=0 


The  orthogonality  relation  for  Q^tx)  is 


n  ( n  +  1 ) (n  +  2 ) 

l  Qk(x,n)QA(x,n)  =  (n  +  ,  _  k)  .(2k  +  1 ) 
x=0  k 


so 


S 
n ) 


I 


k-0 


2 


2k  +  1 


( n  +  1  -  k>k 
<n  +  2)k 


n  2 

Sinr^  (n  +  1  -  k)  /(n  +  2)  <  1,  S  <  \  c  /(2k  +1),  which  proves  one  of  the 

K  k  n  ,  _  K 

k=0 

required  inequalities. 


-6- 


m 


Conversely,  if  Sn  =  0(1),  then 


V  C*  V  ck 

k=0  2k  +  1  -  ^  (2k  +  1) 


( n  +  1  -  k) , 


(n  +2), 


(n  +  2), 


(n  +  1  -  k), 


( n  +  2  )  n  c, 

<  - .  v  — K 


(n  +  1  -  n)  k=0 
n 


2 

k _ 

2k  +  1 


(n  +  1  -  k), 


(n  +  2), 


(1 


2  ‘ 
n  +  1 

< - _s2  =0(1)  . 

n  +  1 


3 .  A  Weighted  Hausdorff  Moment  Problem 

Extensions  of  Legendre  polynomials  and  the  discrete  Tchebychef  polynomials 
exist,  so  it  is  natural  to  see  if  they  can  be  used  to  obtain  an  extension  of 
Hausdorff' s  theorem.  To  this  end,  set 


(3.1)  Wv  =  /  f(x)xV  a(1  -  x)^dx,  v  =  0,1,2,...  (a;P  >  -1)  . 

0 

ci  8 

Polynomials  orthoaonal  with  respect  to  x  (1  -  x)  on  [0,1]  are  known.  Thev  are 
called  Jacobi  polynomials  and  are  given  by 


(3.2) 


P(a'6><1  -  2x) 
n 


(a  +  1  ) 


f-n. 


n+a+8+ 1 
a+1 


Set 


f(x;  ~  )  c  P(a,f>)(1  -  2x ) 

o  v  v 


-7- 


where  c  is  determined  by 


(3.3) 


°v  ~  h(a.0>  0 


~7  /  f(x)P(a,6>{1  -  2x)xa(1  -  x)^dx 


and 


(3.4) 


r 1  -,(a,8),,  .  8>(,  0  v6.  .  , (a,8) 

I  P  (1  -  2x)P,,  (1  -  2x)x  (1  -  x)  dx  =  6,  h 

^  n  k  kn  n 


with 

h«x,B) 

n 

As  in  the  last  section 


_ F ( n  +  a  +  1  )T(n  +  8  +  1 ) _ 

(2n  +  a  +  6  +  1 )T(n  +  a  +  8  +  1 )  •  n! 


.n-v 
A  u 


1 

r  ,,  .  v+a,  „  .n-v+6^ 

J  f(x)x  (1  -  x)  dx 

0 


n  1 

V  r  o  X  v+a,,  .n-v+B, 

*  l  c  ]  P  (1  -  2x)x  (1  -  x)  dx 

k*0  0 

r ( v+a+1  )r(n-v+B+1)  ?  *  1  *k  _  r-k,  k+a+8+1 ,  v+a+1 

r  ( n+a+B+2 )  k^Q  k!  %.  ’  3F2^  a+1,  n+a+B+2 


Using  (2.5)  gives 


I"(n  +  a  +  8  +  2  )An_vy 

_  V 

r<v  +  a  +  i>r(n-v  +  e  +  i) 


(3.5) 


n 

=  l 

k=0 


(a  +  1)] 
ki 


(n  +  1 


k). 


k  (n+a+B  +  2), 


3F2' 


-k ,  k+a+$+1 
a+l ,  -n 


The  aeneral  discrete  Tchebychef  polynomials  {6]  (or  to  use  their  common  namo,  tl’e 
Hahn  polynomials)  are  given  by 


-9- 


Q  (  x ;  a  ,  6  ,  n )  =  F  '  ~k  '  ^ 1+5  +  1  '  _x;  1  ,  (k,x  =  0,1, ...,n>  . 

x  3  2  ci+ 1  ,  -n 


(3.6) 

Their  orthogonality  relation  i- 


r  • x+a  - N-x+ 6 

l  Q.  (  x;a  ,S  ,n)Q  .  (  x;a  ,3  ,  n  )  v 

x=0  K 


N-x 


(3.7) 


(a  +  S  +  2) 


=  6 


k.'(n  +  a  +  3+  2),(3  +  1)  (a 
k  k 


+  1 


jk  n!(n  +  1  -  k)v(o  +  1)^(cc  +  °  +  1)v(2k  +  a  +  3  1) 


k 


0  <  j ,  k  <  n  . 


Square  (3.5),  multiply  by  and  sum.  After  simplification, 

resulting  identity  is 


the 


v  iAn_vii  i2rn) - iiel. 

un  v  ^  v'  r  (u  +  a  +  1  )  F(n  -  v  +  1 )  ~k'  k 

u=o  k=n 


r  ( n  +  a  +  B  +  2) 


n  -  ,  . .  ( n  +  1  -  k ) , 

c  2,  ( a ,  2  )  k 

C  •  - - - 

( n  +  ot  +  £  +  2 ) , 


The  Riesz-Fisher  theorem  for  Jacobi  series  is 


f  (f(x))2xa(1  -  x)Bdx  = 


k=0 


„2„(a,a) 

~k  k 


so  an  argument  similar  to  the  one  in  §2  gives  the  following: 
Theorem  1 .  Define  b^_  (3.1).  Then  for  a,r  >  -1, 

1 

/  |f(x)|2xa(1  -  x/  vix  <  •» 

0 


if  and  only  if 


n 

v 


V=0 


|A 


n-v 

U 


v 


n '  _ r  ( n  +  a  *  ::  *  2  ) _ 

v;  f  (  v  +  a  +  1 ) f ( n  -  v  +  n  +  1) 


0(1)  . 


This  can  hie  rephrased  as 


B  ,  2 


/  |  f  (  x)  x  ( 1  -  x)  | 


dx 


a,, 

x  (1 


-  x)‘ 


if  and  only  if 


v  + 


n  + 


_1_ 

2 

1 


when  a ,  8  >  -1 ■ 


=  0(1) 


-10- 


REFERENCES 


(1]  R.  E.  Bellman,  R.  E.  Kalaba  and  J.  A.  Lockett,  "Numerical  Inversion  of  the 
Laplace  Transform",  Amer.  Elsevier  Publ .  Co.,  Inc.,  New  York,  1966. 

[2]  Felix  Hausdorff,  Moment  probleme  fur  ein  endliches  Intervall,  Math.  Z.  16 
(1923),  220-248. 

[3J  I.  J.  Schoenberg,  Remarks  concerning  a  numerical  inversion  of  the  Laplace 

transform  due  to  Bellman,  Kalaba  and  Lockett,  J.  Math.  Anal.  Appl.  43  (1973), 
823-828. 

[4]  I.  J.  Schoenberg,  On  finite  and  infinite  completely  monotonic  sequences.  Bull. 
Amer.  Math.  Soc.,  Feb.  1932,  72-76. 

[5]  G.  Szego,  "Orthogonal  Polynomials",  Fourth  edition,  Amer.  Math.  Soc., 
Providence,  RI ,  1975. 

[6]  P.  L.  Tchebychef,  Sur  1 ' interpolation  de  valeurs  eauidistants,  Oeuvees  de  P.  L. 
Tchebychef,  Vol.  II,  reprinted  by  Chelsea,  New  York,  1961,  217-242. 


PA/I JS/AS/scr 


SECURITY  Cl  ASGKU  AT:  *n  r 


REPORT  DOCUMENT/*!  ION  PAGO 


1.  REPORT  r«UMt‘ILD 

2130  ^ 


M  \ I .• 
m  m. h.:  r  < 


4.  TITLE  (nnj  Subtitle) 

HAUSDORFF’S  ' 

MOMENT  PROBLEM  AND  EXPANSIONS 

IN  LEGENDRE  , 

.  ‘ VNOMIALS 

7.  AUThORCs) 

R.  Askey,  I. 

J.  Schoenberg,  and  A.  Sharma 

>j  type  or  ."-i  po at  r.  fit  :  r 

Summary  Report  -  r.o 
rn;  ortina  r.f  :;o'i 

C  PERFORMING  O«0.  «Lf  CUT  i 


>.  performing  orgav:ation  name  anti  address 

Mathematics  Research  Center,.  University  of 
610  Walnut  Street  Wisconsin 

Madison.  Wisconsin  33706 


1 1.  CONTROLLING  OFFICE  name  and  address 

U.  S.  Army  Research  Office 

12.  REPORT  DATE 

October  1980 

P.O.  Box  12211 

Research  Triangle  Parle.  North  Carolina  277  09 

13.  NUMBER  OF  PAGES 

11 

14.  MONITORING  -GENCY  NAMES  AODPESSfJ/  different  from  Controlling  Oilicc) 

15.  SECURITY  CLASS,  (ct  t!  us  re;  or:. 

UNCLASSIFIED 

15*.  DECLASSIFICATION  D OWN GR * C  1  *  C* 
SCHEDULE 

16-  DISTRIBUTION  STATEMENT  (of  this  Report) 

Approved  for  public  release;  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  in  BfocJc  20,  if  different  fro. 

m  Report ; 

18.  SUPPLEMENT  ARY  NOTE  S 

19.  KEY  WORDS  (Continue  on  reverse  side  if  necessary  and  Identify  by  block  r-jtr.ber) 

Moment  problem 

Orthogonal  polynomials 

20.  ABSTRACT  (Continue  on  -eveese  side  If  necessary’  and  identify  by  block  number} 

\A  new  proof  is  given  for  Hausdorff's  condition  on  a  sot  of  moments  which 

determines  when  the  function  generating  these  moments  is  in  LG-  The  proof 

uses  Legendre  polynomials  and  their  discrete  extensions  found  bv  Tehobvehef. 

•  •  2  4 

Then  an  extension  is  given  to  a  weiqhted  Lj3  space  using  Jacobi  polynomials 

and  their  discrete  extensions. v 

DD  ,  jaTI-3  1473  EO.HON  OF  1  NOV  6*  IS  OBSOLETE  UNCIASS1  HKD 

SECURITpCl  ASSirTc  A1  ION  Of  THIS  PAGE  ft*!.i-n  I  Mir 


